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Knot-controlled ejection of a polymer from a virus capsid 
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We present a numerical study of the effect of knotting on the ejection of flexible and semiflexible 
polymers from a spherical, virus-like capsid. The polymer ejection rate is primarily controlled by 
the knot, which moves to the hole in the capsid and then acts as a ratchet. Polymers with more 
complex knots eject more slowly and, for large knots, the knot type, and not the flexibility of the 
polymer, determines the rate of ejection. We discuss the relation of our results to the ejection of 
DNA from viral capsids and conjecture that this process has the biological advantage of unknotting 
the DNA before it enters a cell. 

PACS numbers: 87.15.-v, 82.35.Lr, 02.10.Kn 



Icosahedral bacteriophages are viruses that infect bac- 
teria. They typically consist of an almost spherical cap- 
sid head with dimensions of several tens of nanometers, 
on the order of the persistence length of DNA, and a 
narrow cylindrical tail with an internal diameter of only 
a few nanometers, through which the phage injects its 
DNA into bacteria [l[ . Remarkably long strands of DNA 
can be packed to almost crystalline densities inside the 
rigid capsid heads. For example, A phage has a genome of 
length 16 microns squeezed into a spherical capsid with a 
diameter of just 58 ran. Internal pressures can thus be on 
the order tens of atmospheres. A number of viruses, such 
as as A and 029 phages, exploit this pressure to force their 
DNA through their tail, and into their bacterial hosts. 

Experiments, using fluorescent staining and light 
scattering [3, 0], have recently investigated DNA ejec- 
tion from viral capsids. These show that the ejection 
rate can be affected by temperature [3, 0|, the presence 
of binding proteins [H , genome length [i| and the concen- 
tration of salt or other ions Further work has shown 
that ejection can be suppressed, so that only a fraction 
of the genome is emitted, by adding polyethylene glycol 
to change the osmotic pressure of the solution surround- 
ing the capsid In addition, pauses, which may be 
correlated with position along the chain, are seen in the 
ejection of certain phages @. A number of the generic 
features of these experiments can be explained by treat- 
ing the DNA as a simple model polymer which is driven 
from the capsid by the energetic and entropic penalty of 
close confinement [1, 0, 0, S G3 • 

In an intriguing set of experiments, Arsuaga et al. [HI, 
Il2j ]. directly extracted DNA from tailless bacteriophage 
mutants and showed that it was highly knotted, sug- 
gesting that this may also be the case inside the cap- 
sid. Moreover, recent simulations have provided further 
evidence for the prevalence of knots on a confined poly- 
mer: Micheletti et al. [H| showed that the probability of 
knotting for a polymer contained in a sphere increases 
with the polymer length and the degree of confinement, 
in agreement with earlier work [ijj]. Although the spec- 
trum of knots seen in experiments was not exactly the 
same as those from simulations of a random polymer, a 
fact which presumably reflects the chiral nature of the 



DNA packing Q3, the fact that DNA is highly knotted 
is remarkable. Knots can prevent the transcription of 
DNA by RNA polymerase, and cells have developed a 
number of active ways to control these entanglements, 
for example through the use of molecular machines like 
topoisomerases [l|. It seems unlikely that viruses would 
use topoisomerases to unknot their DNA, so instead a 
different mechanism must be at work. 

Motivated by these experiments, our aim in this letter 
is to study the effect of knots on the ejection of a viral 
genome from a phage capsid. The head-tail connector 
has a channel of diameter 2.3 nm in A [l5| and 1.7 nm in 
4>29 [l6[ , compared with an interstrand spacing of about 
2.5 nm for packed DNA [l7|, HI]. Considering that for 
a knot to pass multiple strands must go through simul- 
taneously, even simple knots in DNA are expected to be 
too large to fit. We assume that instead the DNA must 
be extruded by reptating through the knots. A biological 
advantage would be that the viral DNA would enter the 
bacteria unknotted. 

Building on the success of previous modeling [9l. IToL fT9|| 
that reproduced generic effects seen in experiments, the 
approach we take is to represent the DNA by a simple 
bead-spring polymer. The polymer is initially confined 
to a spherical capsid and is coupled to a coarse-grained 
solvent. It is allowed to eject through a small hole, and 
the driving force is the pressure of the packed DNA in the 
capsid. We find that the knots control the rate at which 
the polymer can leave the capsid, with slower ejection 
rates observed for more complex knots. The knot acts 
as a ratchet, with the polymer being ejected as the knot 
reptates along it. 

The coarse-grained polymer we consider comprises 
beads connected linearly by springs. The beads interact 
through the potential 
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FIG. 1: Schematic of initial configurations: Left, a spooled 
chain, right, a random configuration is achieved by packing 
the chain through a hole at the bottom of the capsid, which 
is closed before equilibration. 

The first term in Eq. fl} is the FENE spring potential, 
the second a repulsive Lennard-Jones term representing 
an excluded volume interaction between the beads and 
the final term is a bending potential which can be used 
to control the persistence length of the polymer. We 
choose e = a = 1, k = 30 and i? = 1-5. We consid- 
ered both flexible chains, k — 0, and semiflexible chains, 
k = 10, which corresponds to a persistence length of 10 
beads, on the order of the diameter of the capsids, as is 
the case for real bacteriophages. The dynamics of the 
beads was simulated by using a velocity Vcrlct molecular 
dynamics algorithm. The viral capsid was modeled as 
a hard spherical shell by applying a force of magnitude 
ksT / {a f A ) to beads when their positions satisfied the 
inequality / < 0.2, where / = 1 — {x 2 + y 2 + z 2 )/R 2 . 
In addition, we added a single hole small enough that 
only one bead at a time could pass through. Two sizes of 
capsid were considered: one with a radius of R = 3.02cr 
and a second with R — 4.36<r, which leads to a volume 
three times larger. The polymer was coupled to a sol- 
vent modeled using a stochastic rotation dynamics algo- 
rithm 20]. This provides a thermostat, which conserves 
momentum, and hence means that hydrodynamic inter- 
actions between polymer beads are included in the sim- 
ulation. The capsid was permeable to the solvent, the 
physical case for phage capsids. 

Two types of initial configuration, shown in Fig. [TJ 
were considered. In both cases the knot was put in by 
hand in a tight configuration near the exit, its type con- 
firmed using the Alexander polynomial. In the first case 
the remaining beads were positioned to form a spool [ljj • 
In the second they were initialized in a random configura- 
tion outside the capsid and the polymer was then packed 
by a motor pulling it through a second hole in the capsid, 
opposite to the first. After the packing was completed, 
the second hole was closed. In both cases the polymer 
was equilibrated (within its local minimum), with the 
first bead held in position just outside the capsid, before 
the start of the ejection. 

Polymers of length 100 beads were used with the 
smaller capsid, and polymers of length 300 and 230 were 
used in the larger capsids for the random packing and the 
spooled packing respectively. All reported simulation re- 
sults are averaged over at least 50 independent runs. 

Fig. [5] shows typical results for the fraction of beads 
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FIG. 2: Fraction of beads remaining in the capsid as a func- 
tion of time for unknotted flexible (•) and semiflexible (o) 
polymers and for flexible polymers with 3i (A), 4i (■) and 6i 
(♦) knots. The knots have a clear effect on the shape of the 
ejection curve, and the rate of ejection is primarily determined 
by knot type. 



left in the capsid as a function of time from the release 
of the polymer. For this figure we used a small capsid 
and a random initial configuration but very similar re- 
sults are found for other packings, initial conditions and 
capsid sizes. We compare results for unknotted chains 
to those with the knots 3i, 4i and 6i. (Here we use 
the notation Ck, where C gives the minimal number of 
crossings in a projection of the knot onto a plane and k 
is a standard way to distinguish between knots with the 
same number of crossings [2lj] .) Fig. [2] shows that there 
is a clear slowing of the rate of ejection when a knot is 
present. Moreover, the rate of ejection depends on the 
type of knot: the more complex the knot, the slower the 
rate of ejection. 

At very early times the rates of ejection are similar and 
relatively high for both knotted and unknotted polymers. 
This corresponds to the knot being tightened and pushed 
to the capsid entrance as any free beads between the knot 
and the exit are ejected. The knot is too large to escape, 
so once it has moved to the hole it is held there by the 
excess pressure inside the capsid. Now the polymer has to 
reptate through the knot before any monomers are free to 
allow further ejection. Essentially the knot is acting as a 
ratchet. If it diffuses a small distance into the capsid, this 
frees a length of polymer between the knot and the capsid 
entrance. The knot is then quickly pushed back towards 
the entrance by the driving force and the free section of 
polymer is ejected. Near the end of the ejection, when 
only about ~ 30 — 40 beads remain in the capsid, the 
ejection speeds up and there is a shoulder in the curves 
in Fig.O This corresponds to the knot becoming undone. 

To confirm the effect of the knot on the dynamics, 
Fig. [3] compares similar ejection curves for a semiflexi- 
ble polymer for length TV = 100 for three different initial 
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FIG. 3: Fraction of beads remaining in the capsid as a func- 
tion of time for semifiexible polymers: (i) unknotted polymer 
(o) (ii) polymer with a 6i knot initially at the mid-point of 
the polymer chain (A) and (iii) a polymer with a 6i knot 
initially near the capsid entrance (□). 



FIG. 4: Average ejection rates for Ci knots: flexible polymers 
in the smaller capsid (•); semifiexible polymers in the smaller 
capsid (A); flexible polymers in the larger capsid (■); semi- 
flexible polymers in the larger capsid ((>)• For more complex 
knots, the rates for flexible and semifiexible polymers are very 
similar. 



conditions (i) unknotted, (ii) with a knot initially at the 
the polymer mid-point (iii) with a knot initially at cap- 
sid exit. For case (ii) the polymer initially ejects with 
the same speed as the unknotted polymer (i), but, once 
the knot has reached the capsid exit, the rate becomes 
similar to case (iii). 

Simulations for longer polymers in a capsid of radius 
4.36cr showed the same generic behavior. Some small 
differences were that the final unknotting occurred with 
a somewhat larger number of beads left in the capsid 
because there was more free space, and that the rates 
sometimes showed a slight decrease with time, most likely 
because the pressure decrease has an effect on knot dy- 
namics. For the N = 300 randomly packed polymer we 
sometimes also observed initial jamming, with large vari- 
ations in ejection time between runs. This occurred be- 
cause the knot was initially trapped in a very tight con- 
figuration. However once ejection did start, the rate was 
independent of the initial jamming time. 

A quantitative comparison of ejection rates is shown 
in Fig. |H Results are presented for both flexible and 
semifiexible polymers for knots with up to 6 crossings for 
N = 100 chains in the small capsid, and for knots with up 
to 12 crossing for N — 230 chains in large capsids. The 
results plotted are for the knots of type C\. Other knots, 
such as 52 and 62, were considered. The results were very 
similar to the C\ knots with the same C. Average rates 
were obtained by measuring the time taken for between 
75% to 35% of the beads to be expelled. This protocol 
was chosen to avoid any early time jamming and the late 
time unravelling of the knot. Small changes in the cut- 
offs gave no qualitative changes in the results. 

Fig. 0] shows the trend of decreasing ejection rate with 
increasing knot complexity. For the unknotted case, 
there are clear differences between the different kinds of 



polymers and capsids, but for more complex knots, the 
rates converge and seem to be almost completely deter- 
mined by the knot type. This suggests that, for tight, 
confined knots, the rate of reptation of a knot may be 
independent of its flexibility. A similar independence of 
knot properties - the minimum len gth -to-radius ratio re- 
quired for a tube forming a knot [22| and the diffusion 
constant [23] - from flexibility has been reported. 

Finally, note that the ejection rates show a weak os- 
cillation as a function of knot size. The plausible expla- 
nation for this is that the knots we consider here with 
even and odd crossings belong to families with different 
topologies. Knots 3i, 5i and 1\ are torus knots and 4 1; 61 
and 81 are even twist knots. It is reasonable that there is 
a weak dependence of the reptation rate along the chain 
on the knot topology. We have observed a similar os- 
cillation of the polymer diffusion coefficients of knots on 
unconfined polymers under tension (24|. 

Our model of a bacteriophage is highly simplified. For 
example, we ignored the tail. We checked that this does 
not have an important effect on the qualitative behav- 
ior by simulating a capsid with a tail attached. Per- 
haps more importantly, to make simulation feasible, our 
model of DNA inevitably ignores much chemical detail 
and our polymer diameter to persistence length ratio is 
much larger than in real DNA. Moreover, viral DNA is 
typically much longer than the polymers we study. How- 
ever, we argue that the generic effects we observe: that 
the ejection rate is dominated by the rate of knot repta- 
tion, and that the knot is unravelled as it emerges from 
the capsid, are basic properties of knotted polymers mov- 
ing through a narrow exit hole, and should therefore be 
robust to the inclusion of these details. Knots in phages 
may initially be loose but we expect that they will be 
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FIG. 5: Fraction of beads remaining in the capsid as a func- 
tion of time for semifiexible polymers with three knots ran- 
domly selected from 3i, 4i, 5i, 52 and 6i (A) and the prime 
knot 12i (o). 



tightened as the DNA is forced through the capsid en- 
trance. 

Experiments suggest that there could be multiple knots 
in the DNA in a viral capsid [ll]> [l^] • To investigate the 
effect of multiple knotting we ran 50 simulations in each 
of which three knots, chosen randomly from 3i, 4i, 5i, 
52 and 6i were placed on the chain, near to the entrance 



of a larger capsid. Fig. [5] compares the averaged ejection 
curve to that for a single 12 1 knot, showing that the 
knots stack up at the capsid entrance and behave like 
a single prime knot with a similar number of crossings. 
In the longer DNA in a viral capsid one might expect 
to see multiple knots which would sequentially collect at 
the exit. Thus the knot length would increase with time, 
which is expected to slow the ejection rate. It will be 
interesting to see if any sign of this effect can be seen 
experimentally. 

To summarize, we have used simulations to determine 
the effect of knotting on the ejection of flexible and semi- 
flexible polymers from a spherical capsid. We find that 
the ejection rate is controlled primarily by the knot, not 
the pressure of confinement. The knot moves to the hole 
in the capsid and then acts as a ratchet. Polymers with 
more complex knots eject more slowly. For tightly con- 
fined knots the flexibility of the polymer is not key in 
determining the rate of ejection. 

Repeating packing experiments [25[ with knotted DNA 
would shed light on the dynamics of knots on polymer 
chains, and it is also of interest to ask whether knots 
can affect the motion of biomolecules as they traverse 
nanopores |26| . 
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